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Abstract 

We study the stability of the 4-site plaquette order in a model recently 
introduced by Katoh and Imada to explain the spin gap observed in CaV40g. 
We show that all relevant types of order (4-site plaquette, dimerized, Neel) 
can be described within the framework of Schwinger boson mean field theory, 
and we make predictions regarding their relative stability. The results are 
compared to exact diagonalization of finite clusters. 
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In a recent paper, Taniguchi et al. [0] reported the presence of a spin gap in CaVziOg. 
The structure of the two dimensional crystal of VO4 pyramids of CaX^Og is shown on 
figure [I]. In the Mott insulating phase, one can consider Vanadium atoms as localized 
spin one half coupled by superexchange interactions, and, as pointed out by Katoh and 
Imada |2j a minimal model consists of spin 1/2 antiferromagnetic Heisenberg model including 
an exchange term J between nearest neighbors Vanadium atoms. These authors have shown 
using Monte Carlo simulations that this model has a spin gap. In order to determine the 
origin of the gap, they introduce a slightly more general Hamiltonian: 

h = jy j s 1 .s 3 + j' J2 (i) 

<ij> <ij>' 

where <ij> (resp. <ij>') means full line (resp. dashed line) links on figure |l| Using first 
order perturbation theory in J'/ J, they argue that the system keep the same type of order 
as for J' = i.e. a product of singlet on each plaquette. 

In fact, introducing two links gives us a new quantum spin model which might have 
interesting phase transitions: For J 3> J' the system is gapped and the ground state is 
the product of the singlet state on each square plaquette while for J' J the origin of 
the spin gap is quite different. The ground state is then the product of singlets on each 
inter-plaquette bonds. The question that naturally arises concerns the way the system is 
passing from one type of order to the other while increasing J' /J. Another question that 
can also be asked is whether Neel long range order can exist for S = 1/2 - the system is 
not frustrated! - and, if this is the case, what is the nature of the transition between this 
ordered state and the previously mentioned dimerized states. 

In this Rapid Communication, our aim is to answer these questions using a large S 
approach. A linear spin wave (LSWT) calculation was done very recently by Ueda at al. [|3] . 
This theory predicts that the system is Neel ordered for a large region of parameters: .25 < 
J' I J < 6. As these authors say, such the bounds seem to be incompatible with second order 
perturbation theory in J'/ J and J / J'. This could be explained by the tendency LSWT to 
favor Neel order. Still, the perturbation theory is not expected to be reliable around J' = J, 
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and the possibility of stabilizing Neel order in that region remains. 

Another method to check the presence of Neel order in the context of a large S approach 
is to use the Schwinger boson mean field theory (SBMFT) |J. One advantage is that it is 
a improvement over LSWT because it includes higher order corrections in a self-consistent 
way. But more importantly, this method is able to described not only the Neel ordered 
phase but also the two types of valence bond order which are present in the limits J — > 
and J' — > 0. This method starts from a representation of the spin algebra in terms of bosonic 
operators: Si = \b\ a d ! aa 'bia' , the size of the spin being fixed by a constraint on the number 
of particles on each site: b\*b^ + foji^ij. = 2S\ Defining operators that are quadratic in terms 
of the bosonic operators by 2B}j = b\+bji + b\<bji and 2A\j = b\^bj, — b\,bjp the Hamiltonian 
can be written: 

H = Y,Jn{--Bt j B ir .-A\ j A ij ) (2) 

(id) 

At the mean-field level, one introduces the following order parameters: < A\a >= 2a%j and 
<Bjj>= 2flij and the Hamiltonian is replaced with: 

Hmf = E J ij (Aj (Bij + 4) " <*ij (Aj + 4) - Pfj + 4) ( 3 ) 

Finally, the local constraint is replaced by a global one and is enforced only on the aver- 
age through the addition to the Hamiltonian of a term (^It^T + ^IjAi — 25^, where 
the chemical potential fi plays the role of a Lagrange parameter. Because there are only 
antiferromagnetic like correlations in this model, we have: < Bij >= for all pairs of spin 

(hj)- 

These equations have two types of solution: In the first one, the dispersion relation is 
vanishing at three points in the Brillouin zone. This means that the bosons are condensed 
with a macroscopic expectation S* for the occupation number of the zero energy modes with 
the consequence that 

lim < Si.Sj >= ±S* 2 (4) 



The second kind of solution has a gap in the excitation spectra. There is then no long range 
order. The mean field equations have two remarkable solutions of this type: one has the 
same energy and the same correlation functions as the plaquette-singlet one, while the other 
as the same physical properties as the dimerized one. 



SBMFT approximation, there is only one solution for a given value of J' / J. For J' / J < 0.6 
we find a gapped solution with an energy per site —.5 J: this solution has the same correlation 
functions as the plaquette-singlet. At the critical value J' /J ~ .6 the gap closes. For 
.6 < J' I J < 2.4 the system exhibits Neel long range order. At J' /J = 2.4 a gap opens 
again and the order corresponds to dimers on the J' bonds. The order parameters < > 
are continuous through both transitions which shows they are second order. Note that this 
result is meaningful and not a mere consequence of the approximation used. For another 
model, using SBMFT, we found a first order transition between LRO and valence bond 
order |2J. 

We have also done exact diagonalisation of small clusters in order to check our results 
and to have an estimation of the gap in the spin excitations. As there are four atoms per 
unit cell, one can only study square clusters having AN atoms. Besides if N is odd, periodic 
boundary conditions are frustrating. So the only tractable clusters have 8 and 16 spins. 
The next cluster (24 sites) dose not have the symmetry of the square lattice. For these two 
clusters, we have calculated the ground state energy, the spin gap A^ and the staggered 
magnetization defined by M* =< (j)\(l/NJ2i e i^i) 2 \4> > where \<p > is the ground state of 
the system and = ±1 depending on the orientation of the corresponding classical spin in 
the Neel state. 

The different quantities can be scaled to their thermodynamic limit ||: 



The phase diagram of the model under investigation is represented on figure |3|. At the 




AT3/2 



(5) 



A N = A + -f 



(6) 
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M n = M; + j^ (7) 

The results are displayed on figures |2||| and |j. Taken literally, the results are inconsistent 
because the magnetization and the gap cannot be non zero at the same time. But given 
the sizes we could look at, very small values of M* or A are not accessible and only sizable 
values should considered as meaningful. With this proviso, we find a qualitative agreement 
between exact diagonalisation and SBMFT. The staggered magnetization we deduced from 
finite size extrapolation is rather large for .6 < J' /J < 2.4 and comparable to the value 
calculated within SBMFT. In that region of parameters, the value of the extrapolated gap 
is small and reaches a minimum of 0.04 J at J'/ J = 1.25. This two facts suggests that the 
SBMFT calculation are reliable, and that the system under investigation sustains Neel long 
range order for J'/Jza 1.25 

To summarize, the plaquette singlet order proposed by Katoh and Imada to explain the 
magnetic properties of CaX^Og can be found in the context of 1/S approach provided one 
goes beyond LSW theory by using SBMFT [JF]). This method enables one to study the 
stability of this type of order with respect to the switching on of a coupling J' between the 
plaquettes. The result is that Neel order is recovered in a sizable region around J' /J ~ 1.25, 
to disappear again when J' / J is large enough in favor of a dimerized order. It would be very 
interesting to see if Monte Carlo calculations of the gap sustain the presence of a gapless 
region. 

At a more quantitative level, let us note that the predictions of SBMFT disagree with 
the results of Katoh and Imada, who reported a gap of about 0.1J for J' = J, while the 
SBMFT predicts the system to be already gapless for J' / J « 0.6. A possible explanation is 
that there is indeed an ordered phase but it is smaller then what SBMFT predicts, and in 
particular it does not include the point J' = J. Note however that, contrary to Katoh and 
Imada's results, very recent Monte Carlo simulations by Troyer et al. || suggest that Neel 
order is present for the case J' = J. More work is needed to check that point. In any case, 
as far as CaX^Og is concerned, this might be of minor importance. As pointed out recently 
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by Sano et al. ||, the plaquette order can be stabilized by introducing some frustration into 
the system. 

We thank the authors of reference ||, and especially H. Kontani, for pointing out to 
us a mistake in the first version of this manuscript. We also thanks Matthias Troyer for 
informing us of his results prior to publication. The numerical calculations reported here 
were possible thanks to computing time made available by IDRIS, Orsay (France). 
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FIGURES 

FIG. 1. (a) Schematic structure of the VO4 pyramids plane, (b) effective Heisenberg model 

FIG. 2. Ground state energy J'/ J for N = 8 (o) and N = 16 (o) clusters and values extrapo- 
lated for N = 00 (•) together with SBMFT energy (full line) 

FIG. 3. Magnetization versus J'/ J for N = 8 (o) and N = 16 (o) clusters and values extrapo- 
lated for N = 00 (•) together with SBMFT energy (full line) 

FIG. 4. Spin gap versus J' / J for N = 8 (o) and N = 16 (o) clusters and values extrapolated 
for N = 00 (•) 
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